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Abstract. In this paper we show that there exists a family of supersonic neutral 
modes for a compressible mixing layer in an unbounded domain. These modes have zero 
wavenumber and frequency with non-zero phase speed. They are analogous to the super- 
sonic neutral modes of the compressible vortex sheet found by Miles. The results 
presented here give a more complete picture of the spectrum of the disturbances in this 
flow. 
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1. Introduction. In recent studies (Jackson and Grosch, 1989, 1990), we considered 
the inviscid spatial stability problem for the compressible mixing layer in an unbounded 
domain. The key to understanding the stability characteristics of this flow is the under- 
standing of the different parameter regions for which various types of instability modes 
can exist. In Figure 1 we show a plot of the sonic phase speeds c ± versus the Mach 
number, where 


c+ 



c_ - Py + 


VPr 

M 


Here, M is the Mach number of the fast stream, Py the ratio of the speed of the slow 
stream to that of the fast stream, and JJj- the ratio of the temperature of the slow stream to 
that of the fast stream. These curves divide the phase speed-Mach number plane into four 
regions. If a disturbance exists with a Mach number and phase speed in region 1, the dis- 
turbance is subsonic at both boundaries, and we classify it as a subsonic mode. In region 
3, the disturbance is supersonic at both boundaries, and we classify it as supersonic- 
supersonic mode. In region 2, the disturbance is subsonic in the fast stream and super- 
sonic in the slow stream, and we classify it as a fast mode. Finally, in region 4, the dis- 
turbance is supersonic in the fast stream and subsonic in the slow stream, and we classify 
it as a slow mode. We found that there is only a single subsonic neutral mode for two 
dimensional waves. Beyond a critical Mach number M s , the Mach number at which the 
phase speed equals that of a sonic wave, the subsonic neutral modes are transformed into 
supersonic neutral modes which are subsonic in one freestream and supersonic in the 
other. In addition, another supersonic neutral mode appears at M* ^ M s , the Mach 
number at which the sonic speeds of the fast and slow streams are equal. This supersonic 
neutral mode has the opposite behavior than the previous. That is, if the continuation of 
the subsonic neutral mode is supersonic in the fast stream and subsonic in the slow 
stream, this new mode is subsonic in the fast stream and supersonic in the slow stream. 

In these studies we did not find any neutral or unstable modes in region 3 which 
were supersonic in both streams. In this paper we now show that there exists a family of 
neutral modes in region 3 having zero wavenumber and frequency with non-zero phase 
speed 1 . The approach taken here follows that of Papageorgiou and Smith (1989) and 
Papageorgiou (1990), who considered long waves in the wake behind a flat plate. We 
show that these modes are analogous to the supersonic neutral modes of the compressible 
vortex sheet found by Miles (1958). The results presented here give a more complete pic- 
ture of the spectrum of the disturbances in a compressible mixing layer. 

2. Formulation and Results. We consider here a compressible mixing layer in an 
unbounded domain with mean velocity U (q) and temperature T (q), where q is the simi- 
larity variable in terms of the downstream distance and the Howarth-Dorodnitzyn variable. 
With suitable normalization, U varies between 1 in the fast stream (q = <*>), and < 1 in 


This possibility was suggested to us by R. Klein. 
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Ihe slow stream (T| = -<*). The temperature T varies between 1 in the fast stream and 
Pf > 0 in the slow stream. The results presented here are independent of the detailed form 
of U and T. 

The mean flow is perturbed by introducing two-dimensional wave disturbances. The 
disturbance equation for the normal velocity perturbation <|> is given by (Lees and Lin, 
1946) 


( U - c)4>'- lf± 
G 


-a 2 (U — c ) <(> = 0, 


( 1 ) 


where primes denote differentiation with respect to the similarity variable rj, and G is 
given by 


G =T[T -M 2 {U -c) 2 ]. 


( 2 ) 


Here, ct is the wavenumber of the disturbance in the downstream direction and c is the 
corresponding phase speed. 

Since we are interested in long wavelength perturbations (a — » 0), the appropriate 
expansions are found to be 


(J) = <)> 0 + a4>! + ot 2 4> 2 + • ■ ■ 

(3) 

c - c 0 + ac[ + a 2 c 2 + 

(4) 


—— = + cc G < + ct 2 G 2 + ■ ■ ■ (5) 

G G o 


where 


Gq = T [T — M 2 (U — Cq) 2 ], 

(6) 

G\ = —1c\M 2 Gq 2 (U -c q )T. 

(7) 


Let us now define the linear operator L : 


LQ¥) = {U - Cq)^' - U >X V = (U - c 0 ) 2 


¥ 

U -c 0 


( 8 ) 


Substituting the above expansions into (1) gives the following three leading-order prob- 
lems: 
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d r| 


(Go 1 L(%)) = 0, 


-j- (G o 1 [L $0- c j ^3) + -J- (G j L «|) 0 )) = 0, 

arj «*! 


-f - (G o" 1 [L (<|> 2 ) - C 2 - c ! <M) + -J- (G J [L (fo) - c J <M) 

dr \ uT| 

+ -f - (G 2 L (<J> 0 )) - (U - c q) <J>q = 0. 
dr \ 

The appropriate solution to (9) that is bounded for large Ti is given by 


(9) 

( 10 ) 


( 11 ) 


^0 ~ Ao(U — Cq), 


( 12 ) 


where Aq is a constant. Based on this solution the following three remarks can be made: 

[1] Since solutions in region 2 of Figure 1 are subsonic in the fast stream and 
supersonic in the slow stream, we see that <J>o is a solution if we take the 
leading-order wave speed to be 

c 0 * 1. 0 3 ) 

[2] Since solutions in region 4 of Figure 1 are subsonic in the slow stream and 
supersonic in the fast stream, we see that <|>o is a solution if we take the 
leading-order wave speed to be 

c 0 = $u ■ 

[3] Since solutions in region 3 of Figure 1 are supersonic in both the fast and 
slow streams, we see that <J> 0 is a possible solution for some leading-order 
wave speed lying in the range 

c_ < Cq < c + , (15) 

where c ± are the phase speeds of the sonic neutral modes. Thus, the actual 
value of c 0 is not determined at this order, and can only be found by using 
the higher orders terms in the a expansion. 

To determine the leading-order wave speed c 0 for region 3, we consider the expansions at 
the higher orders. The solutions to (10) and (11) are given by 


Co 

((/ - c 0 ) 2 


d ti. 


<J>! =-c 1 A 0 + A X (JJ - c 0 )J 


(16) 
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. <j»l 

^2 = -c 2 A 0 + c x (U - c Q ) J — -dT\ 

(U - c Q y 

+ A 0 (U -c 0 )/ -- G ° 2 f f(t/ -C 0 ) 2 ds ] dT\. (17) 

((7 - c 0 r L J 

We note here that both <jq and 4»2 could contain another term that is proportional to <J>o» 
which can be ignored by a suitable renormalization of i4g. 

Inspection suggests that the solution diverges as q — » ± °°. Thus, these expansions 
will break down when | T| j = O (a -1 ). There are two outer regions, one in the fast stream 
and one in the slow stream, which will be considered separately. For the outer region in 
the fast stream, the appropriate scalings are 

q = cr 1 z (18) 

for the independent variable, and 


- 4*0 + a ^i + ‘ ‘ ' (19) 

for the dependent variable. Substituting (18) and (19) into our original equation (1) yields, 
to leading order, the equation 


4>o + G$ % = 0, G 0 + = M 2 (l-c 0 ) 2 - 1. 


The solution for outgoing waves is found to be 


= B 0 e 


-i VG o z 


The next order equation is given by 


<f>i + Gq 4>! =gi 4>o, 


8i = 2c 1 A/ z (l - c 0 ), 


whose solution for outgoing waves is 


2vGq 


( 20 ) 


( 21 ) 


( 22 ) 


( 23 ) 


Expanding the solution as z — > 0, and matching with the inner solution (3) as q —> 
yields the following relationships 
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^o = ^oO — c o)> (24) 

B j = — c i A Q, (25) 

A i = i A o ( 1 — c q) 2 (G o ) 1/2 , (26) 

A 2 - i ! (2 - A/ 2 (1 — c 0 ) 2 M 0 - (27) 

(G 0 + ) 3/2 


For the outer region in the slow stream, the appropriate scalings are 

4*5 = ‘ . (28) 

and z is defined as before. Substituting into our original equation (1) yields, to leading 
order, the equation 


$o + G 0 $o _ 0, Go - Pr [M 2 (P(/ - c 0 ) 2 - Py]. 


The solution for outgoing waves is found to be 


<t>o — D o e 


Vg 0 z 


(29) 


(30) 


The next order equation is given by 


<t>i + G 0 <j>i = g\ <t>o, £i - 2c\M 2 $j (Py - c 0 ), 
whose solution for outgoing waves is 


4>l =B>ie 


-i Vg7 z 


+ / 


gi~ D o ■ .-iVST. 

2^ 


(31) 


(32) 


Expanding the solution as z -» 0, and matching with the inner solution (3) as rj -» 
yields the following relationships 


7^0 - Aq($u - Cq), (33) 

£>l = -C\Aq, (34) 


/4 1 - i A o (P(/ — c 0 ) 2 (Gq ) 1/2 » 


(35) 
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A 2 = i 


P u ~ Cp 
(Gq ) 3/2 


Pi M 2 JV -M 2 (P v -C 0 ) 2 )A q . 


(36) 


The two conditions (26) and (35) on A j yields the following eigenvalue equation for 
c 0 : 


P t [W 2 (P u ~ c 0 ) 2 ~ P /■ 1 0 - c 0 ) 4 = [M 2 (l - c 0 ) 2 - 1 ] (P u- c 0 ) 4 . (37) 


This equation is identical to (5.3a) of Miles (1958) if we re-express his result in our nota- 
tion. Miles showed (in our notation), that for region 3: 

fl] A single root of (37) exists for 


M £ A/. 


i + VP/ - 

T^pT’ 


with phase speed 


c o - 


Pt/ + VP t 
1 + VPr 


(38) 


(39) 


which we classify as a constant speed supersonic-supersonic neutral mode. 
Note that this solution is independent of Mach number, and corresponds to the 
phase speed at which the two sonic speeds c ± are equal. 

[2] A double root first appears at 


M = 


(1 + (3/ /3 ) 3/2 
1 -P u ’ 


with phase speed 


Co = 


Pt/ + P / /3 
1 + P r ' 3 ' 


[3] There are three distinct real roots for 


M > 


(i + 

1 - Pt/ 


(40) 


(41) 


(42) 


One of the roots is given by (39), while the other two roots must be found 
numerically. For the special case of Py = 1, these roots are given by 
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Cn = 


1 + P U 


2 M 


i — * „ 1 1/2 

M 2 (l - p l/ ) 2 + 4-4VM 2 (l -p t/ ) 2 + 1 . 


(43) 


The root which corresponds to the (+) sign we classify as a fast supersonic- 
supersonic neutral mode, while that which corresponds to the (-) sign we clas- 
sify as a slow supersonic-supersonic neutral mode. 

In Figure 2 we display the real roots of (37) which lie in region 3 for P^ = 0 and 
various values of P^. The roots are the neutral phase speeds. For completeness, we also 
display the neutral phase speeds which lie in the other regions (Jackson and Grosch, 
1989). The classification scheme is given in the figure caption. 

The two conditions (27) and (36) on A 2 yield the result that = 0. Thus, we have 
A 2 = B ! = D j = 0, and the final solution in region 3 is now given by 


— A o (U — c o) 


1 + i a 


(1 - c 0 ) 2 , G 0 
— 1=~\ T dr \ + 

Vg 0 + j W ~ c o) 


a 2 f — - [ J (U - c 0 ) 2 ds ] dr\ + 

1 (U - c 0 ) 2 1 J J 


(44) 


3. Conclusions. We have shown that the compressible mixing layer has a family of 
supersonic neutral modes in region 3 of Figure 1, which have phase speeds which are 
identical to those of the compressible vortex sheet (Miles, 1958). For the mixing layer the 
wavenumbers of these modes are zero, while those of the vortex sheet are arbitrary. One 
might expect this since the vortex sheet has no natural length scale, while the mixing layer 
has a non-zero thickness. We note that these results can be extended to three dimensional 
disturbances by rescaling the Mach number by cos0, where 0 is the angle of propagation. 
Finally, Artola and Majda (1987, 1989) showed that the nonlinear interaction of the three 
supersonic neutral modes of the compressible vortex sheet can lead to instabilities of the 
flow. We surmise that similar behavior might be exhibited by the interaction of the 
supersonic-supersonic neutral modes presented in this paper. We are currently investigat- 
ing this numerically. 

Acknowledgements. The authors are indebted to R. Klein for suggesting that a 
supersonic-supersonic neutral mode might exist in the mixing layer. 
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Figure 1. Plots of the sonic speeds c ± versus Mach number showing the four regions 
in which different types of disturbances can exist. 



- 10 - 



Figure 2. Plot of the neutral phase speeds as a function of Mach number for Pt/ = 0 
and (a) Pj = 2, (b) P r = 1, and (c) p y - = 0.5. The neutral mode classification is: (1) 
subsonic, (2) fast (supersonic), (3) slow (supersonic), (4) = 1, (5) c ^ = 0, (6) con- 

stant speed supersonic-supersonic, (7) fast supersonic-supersonic, and (8) slow 
supersonic-supersonic. The sonic curves are shown as dashed. 
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